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We demonstrate tunable spin-spin couplings between trapped atomic ions, mediated by laser 
forces on multiple transverse collective modes of motion. A axO-^-type Ising interaction is realized 
between quantum bits stored in the ground hyperfine clock states of ^'^^Yb^ions. We demonstrate 
entangling gates and tailor the spin-spin couplings with two and three trapped ions. The use of 
closely-spaced transverse modes provides a new class of interactions relevant to quantum computing 
and simulation with large collections of ions in a single crystal. 
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Trapped atomic ion systems are a leading medium for 
the generation of quantum entanglement in applications 
ranging from quantum computing and communication [2| 
to quantum simulations @] . There have been several re- 
cent demonstrations of controlled entanglement of few- 
ion systems [3j, most involving quantum bit (qubit) state- 
dependent optical forces 0, [1] . Current effort is devoted 
to the scaling to much larger numbers of trapped ion 
qubits. One promising approach entangles small num- 
bers of ions through a coupling to a single mode of col- 
lective motion in a single zone, and shuttles individual 
ions between different trapping zones to propagate the 
entanglement 0, Q ■ It may be difficult to entangle more 
than a few ions through a single mode, owing to poten- 
tial decoherence or uncontrolled coupling with the many 
spectator motional modes 0- However, it may be possi- 
ble to entangle large numbers of ions through all modes 
of collective motion, where scalability relies on the high 
density of motional modes and the relative insensitivity 
to any particular mode [1, [l^l ■ 

In this paper, we demonstrate the tunable interaction 
between two and three trapped ^^^Yb^ions through the 
simultaneous coupling to multiple closely-spaced trans- 
verse modes of motion 



111 Il2l 1131 1 . In addition to their 



high density, transverse modes can be less sensitive to de- 
coherence compared with conventional axial modes 11 1 . 
We apply an optical spin-dependent force on two ions 
at a frequency that bisects the two modes of transverse 
motion in 1-dimension, resulting in a M0lmer-S0rensen 
fajCaj-type synchronous quantum gate between the ions, 
and we demonstrate Bell state entanglement fidelities of 
greater than 95.9% with a gate time of 37.5^sec. When 
the force operates at a frequency far-detuned from the 
two modes, the motion can be adiabatically eliminated 
from the Hamiltonian, resulting in a nearly pure spin- 
spin coupling. For three ions, we create a GHZ entangled 
state based on a specific synchronous quantum gate driv- 
ing all modes. We finally measure two- and three-spin 



Ising dynamics and extract the effective spin-spin cou- 
plings while varying the detuning of the spin-dependent 
force. For larger numbers of ions, this type of control 
of the form and range of the interaction can be applied 
to the simulation of complex models of quantum mag- 



netism. 



We entangle multiple trapped ^^^Yb^ions, represented 
by the ^6*1/2 hyperfine "clock" ground states in each ion, 
abbreviated by the effective spin— 1/2 states |t)and ||), 
respectively, having frequency splitting wq/Stt — 12.643 
GHz [l3|. We uniformly address the ions with two Ra- 
man laser beams in a geometry where the wavevector dif- 
ference Sk of the beams points along the transverse (x) 
direction of trapped ion motion. Each beam has multi- 
ple spectral components that simultaneously drive blue 
and red motional sideband transitions, following the orig- 
inal scheme of M0lmer and S0rensen [1] , but because the 
transverse modes of the ions are closely spaced, they all 
contribute to the interaction. 

In general, when noncopropagating laser beams have 
bichromatic beatnotcs at frequencies luo±ii, this can give 
rise to a spin-dependent force at frequency /j, [3] ■ Under 
the rotating wave approximation {loq 3> /i) and within 
the Lamb-Dicke limit where Sk{xi) <^ 1, with Xi the 
position operator of the ith ion, the resulting interac- 



tion Hamiltonian is H{t) = fli{6k ■ Xi)ax \iTi^t 11|. 



Here, ax is the Pauli spin flip operator on ion i with Rabi 
frequency and Sk-Xi = J^m 'yj^mCame""^'"* + at„e'"™*) 
in terms of the normal mode phonon operators a,„ and 
al^ at frequency ujm- The Lamb-Dicke parameters rji^m = 
bi,mSk^yh/2Mu!m include the normal mode transforma- 
tion matrix bi^m of the ith ion with the mth normal mode 
[1^, where M is the mass of a single ion. 

The evolution operator under this Hamiltonian can be 
written as j9| 
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where 4(r) = Em["i,m(T)aJ„ - "*,m('^)am]- The first 
term in Eq. ([T]) represents spin-dependent displacements 
of the mth motional modes through phase space by an 
amount 

ai,m(T) = ^^^^^^^^^[^-e'"™^(/xcos^r-iw„iSin/ir)]. (2) 

The second term in Eq. ([1]) is a spin-spin interaction 
between ions i and j with coupling 
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We are interested in two regimes where multiple modes 
of motion contribute to the spin-spin coupling, taking 
evolution times t to be much longer than the ion osc- 
cillation periods (wmT ^ 1) [2^. In the "fast" regime, 
the optical beatnote detuning \x is close to one or more 
normal modes and the spins become entangled with the 
motion through the spin-dependent displacements. At 
certain times of the evolution however, ai^m(T) k, for 
all modes m and the motion nearly factorizes, which is 
useful for synchronous entangling quantum logic gates 
between the spins [20]. For the special case of iV = 2 
ions, both modes in a given direction simultaneously fac- 
tor when the detuning is set exactly half way between 
the modes, or at other discrete detunings, where both 
modes contribute to the effective spin-spin coupling. In 
the "slow" regime, the optical beatnote frequency is far 
from each normal mode compared to that mode's side- 
band Rabi frequency — Wmj ^ r]i,„Mi)- In this case, 
the phonons are only virtually excited as the displace- 
ments become negligible (a^^m «C 1), and the result is 
a nearly pure Ising Hamiltonian from the last (secular) 



term of Eq. ([3]): iJg// = ^Z^i , jcri' ^cri''\ where 



(4) 



In our experiment, we measure the spin dynamics in both 
regimes for multiple ions, showing the effect of multiple 
transverse modes on the interaction. 

We trap -'^^^Yb'^ions in a three-layer linear rf-Paul trap 
2ll | , with the nearest trap electrode a transverse distance 
of 100/zm from the ion axis. An rf potential of amplitude 
~ 200 V at 39 MHz is applied to the middle layer, and the 
top and bottom layers are split into axial segments, where 
static potentials of order 10 — 100 V are applied between 
the end and center segments for axial confinement. The 
ions form a linear chain along the trap z— axis, and the 



axial center-of-mass (CM) trap frequency is controlled 
in the range uj^j^-K = 0.6 — 1.7 MHz so that the ions 
are separated by 2-5 [ivt\. In the transverse x direction, 
the highest frequency normal mode of motion is the CM 
mode, set to LOijl-K ^ 4 MHz. The second-highest is 
the "tilt" mode at frequency L02 = — e^, where the 

trap anisotropy e — cOz/oJi <JC 1 scales the spacing of the 
transverse modes. For three ions, the lowest "zigzag" 
mode occurs at frequency CJ3 — cjiVl — 2.4e^. The other 
transverse (y) modes are sufficiently far away from the x 
modes (frequencies below 3 MHz), and the coupling to 
these modes is suppressed by more than a factor of 10 
by rotating the principal y-axis of motion to be nearly 
perpendicular to the laser beams (2.2]. 

We direct two Raman laser beams onto the ions to 
drive spin-dependent forces, with their wavevector dif- 
ference aligned along the transverse a;-axis of ion motion 
(Sk = k\/2, where k = 27r/A). The Raman beams are set 
to a wavelength of A = 369.75 nm, detuned ~ 500 GHz to 
the red of the ^^^Yb"*" "^81/2 P1/2 transition, and have 
a power of up to ^ 8 mW in each beam. The Raman 
beams are phase modulated at a frequency nearly half 
of the ^^^Yb~''hyperfine splitting with a 6.32 GHz reso- 
nant electro-optic modulator, and each of the two Raman 
beams have independent acousto-optic frequency shifters 
to select appropriate optical beatnotes to drive Raman 
transitions in the ions [l^. The Raman beams are fo- 
cussed to a waist of ~ 30/im, and when their beatnote 
is adjusted to drive the carrier transition at the hyper- 
fine transition ojq , we observe a carrier Rabi frequency on 
each ion of up to fl /27r « 300 kHz with inhomogeneities 
between the ions measured to be under 1%. 

We first initialize the ^^-'^ Yb^ions with 3 ms of Doppler 
laser cooling on the ^5*1/2 — ^ P1/2 transition at 369.53 
nm. We then Raman sideband cool all m modes of 
transverse motion along x to a mean vibrational indices 
of n„i < 1 in ^ 0.5 ms, well within the Lamb-Dicke 
limit. Next, the ions are each initialized to the )J,)state 
through optical pumping [l3]. We then apply the op- 
tical spin-dependent force on the ions for a duration t 
by impressing the bichromatic beatnotes at ± M on 
the Raman laser beams. Afterward, the spin states are 
measured by directing resonant laser radiation on the 
^8i/2iF = 1) — ^ Pi/2{F' = 0) transition, following stan- 



dard state-dependent fiuorescence techniques [17|. We 



repeat the same experiment 500 times and average the 
results. With a single ^^^Yb'^ion, the effective detection 
efficiency of the spin state is 97% using a detection time 
of 0.8 ms. For multiple ions, we measure the probability 
Pn of having n ions in the based on the histograms of 



the collected fluorescence [23] 



For entanglement of = 2 ions through multiple 
modes, we set w^/27r = 0.616 MHz, uji/2tt = 3.5838 
MHz, and W2/27r = 3.5305 MHz. Fig. ^a) shows the 
evolution of the spin populations Pq and P2 when the 
Raman coupling is applied with an optical beatnote de- 
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FIG. 1: (a) Time evolution of the spin populations P„ for 
two ions subject to a laser force tuned to bisect the transverse 
CM and tilt sidebands (detuned by 26.7 kHz from each mode) . 
The ions ideally evolve from |J,J,) to the entangled state |J,J,) + 
lit) after a gate time of Tg = 37.5 (dashed line), (b) 
Time evolution of the spin populations for three ions subject 
to a laser force detuned 9.4 kHz blue of the transverse CM 
mode. The ions ideally evolve from to the entangled 

state liii) + lilt) + Itit) + |TTi> after a gate time of r = 
53.4 /is (dashed line), (c) Measured parity of two (circles) 
and three (squares) ions after the respective gate is applied 
followed by a subsequent global 7r/2 rotation of the spins with 
phase (j). (For the three-ion case, an fixed global 7r/2 rotation 
is added before this to create the GHZ state.) The contrast 
of the sinusoidal oscillations with N(j> give the A*'— particle 
coherence of the state. The populations at Tg in (a) and 
(b) along with the parity contrast in (c) imply entanglement 
fidelities of F2 > 95.9% for = 2 and F3 > 84% for iV = 3 
ions. The lines in (a) and (b) are fits to Eq. ([TJ with no free 
parameters, and in (c) are sinusoidal fits to extract contrast. 



tuning /i/27r = {uji + lo-2)I^-'k = 3.5571 MHz that bisects 
the CM and tilt mode sidebands, 26.7 kHz away from 
each sideband. We adjust the Raman laser power so 
that after the gate time Tg — 2tt/{uji — fi) — 37.5 fis, 
the spin-spin coupling from Eq. ([3]) is xi,2 — 7r/4, 
or ryfi = (wi — yLt)/2-\/2, where the magnitudes of the 
Lamb-Dicke parameters for either ion and either mode 
are all approximately rj = 0.049. At time Tg, we measure 
^0 + Pn = 97.8(1.6)%. To characterize the coherence of 



the entangled state, we apply the same gate followed by 
a global Raman carrier 7r/2-pulse with phase </> relative 
to the gate, and measure the parity J2n(~^)^-Pn of 
resulting state, where are the new populations follow- 
ing the rotation. The contast of the parity oscillating as 
sinTVf/) [circles in Fig. [UJc)] of 94(1)% for TV = 2 indicates 
the fidelity of the two-ion entangled state F2 > 96(2)% 
[2^. After 15 consecutive gate operations (562 /is), the 
entanglement fidelity degrades to 74(4)%, implying a 
single-gate fidelity of 98% exclusive of state preparation 
and detection errors. For N = 3 ions, we set the trap fre- 
quencies to uj,j2Tr = 1.484 MHz, tJi/27r = 3.952 MHz, 
a)2/27r = 3.663 MHz, and ^3/271 = 3.215 MHz, and cre- 
ate a GHZ entangled state by tuning the Raman beat- 
note blue of the CM sideband by (^ — tJi)/27r — 9.4 kHz. 
While the other modes do not contribute nearly as much 
as the CM, they must be considered for the appropriate 
gate time and Raman power for the operation. Here, at 
time Tg = 27r/(u;i — /x) = 53.4 /is, we expect Xi,j = t^/^ 
for all spin pairs when the Raman power is set such that 
in w (tji - /i)/1.95, where rj' = 0.0380 is the 3-ion 
CM Lamb-Dicke parameter. Fig. [IJb) shows the time 
evolution of the state, and after the gate time Tg a sim- 
ple 7r/2 rotation of all the spins ideally creates a GHZ 
state ITTT) + liii)- By applying a subsequent carrier 7r/2 
pulse with variable phase as above and measuring the 
constrast of the parity oscillations and populations [Fig 
[Ijc)], we arrive at a three-ion GHZ entanglement fidelity 
of 84(3)% [Po + Pat = 91(2)%, parity contrast = 77(1)%]. 
The main sources of gate infidelity and decoherence are 
classical laser power fluctuations and spontaneous emis- 
sion from the Raman beams. Additional terms not con- 
sidered in the model Hamiltonian, including slight dif- 
ferential AC-Stark shifts, off-resonant transitions on the 
carrier, and higher order Lamb-Dicke effects are expected 
to contribute less than 2% to fidelity errors. 

Next we set the bichromatic detuning ^ to be far from 
all motional modes so that the evolution is dominated by 
pure spin-spin couplings. In this case, we measure Po(r) 
as well as the average number of ions in state |t)P(T) = 
^^nP„(r) after various durations r of the bichromatic 
force. We can then extract the spectrum of couplings Ji j- 
by observing the secular Fourier components in the time 
evolution [Eq. H])]. 

For N=:2, the observed spin evolution is a simple sinu- 
soid shown in Fig. [2Ia), accompanied by low-contrast 
fast undulations from the residual motional couplings 
[first three terms of Eq. We plot the measured cou- 
plings Ji^2 as the detuning ^ is varied in Fig. [2^c), along 
with that expected from Eq. ^ with no adjustable pa- 
rameters. We determine the sign of the coupling by pre- 
ceding the spin-dependent force with a carrier 7r/2 pulse 
and measuring the phase of the spin evolution. 

For = 3 ions, there is competition between the near- 
est neighbor couplings Ji,2 = J2,3 and the next-nearest- 
neighbor coupling Ji,3, as seen in the time evolution of 
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FIG. 2: (a) Time evolution of the average number of ions 
in |T)for N — 2 ions under the bichromatic force in the far- 
detuned limit, showing the secular oscillation of the Ising spin- 
spin coupling at 4.8 kHz accompanied by small fast modula- 
tions from motional excitation at ~ 56 kHz. (b) Time evo- 
lution of the probability Po of all spins |J,), with N = 3 ions 
under the bichromatic force in the far-detuned limit. Here, 
the fast motional oscillations are negligible, but the two cou- 
plings Ji,2 = ^2,3 and Ji,3 are clearly visible. Solid line is a fit 
to the secular terms in Eq. ^ with an empirical exponential 
decay, (c-d) Measured couplings Jij for two (c) and three (d) 
ions (points) as a function of detuning overlaid with theory 
(lines) from Eq. ^ with no free fit parameters (the overall 
sign in (d) is set from theory). The normal modes positions 
are indicated by the dashed lines. 



Pq in Fig. m^b). By scanning the detuning /i, again we 
see from Fig. [^d) that the form of the couplings can 
be controUed. For example, in regimes where the next- 
nearest-neighbors are anti-ferromagnetic (Ji,3 > 0), the 
system is expected to exhibit magnetic frustration. In 
the limit of very large detuning — w^j 3> wi — Ci.'s), 
we expect a dipolar decay of the couplings Ji^3 = 0^1,2/8 
from Eq. ^ for the equally-spaced ions. 

In summary, the experiment demonstrates how to tai- 
lor spin-spin Hamiltonians in trapped ion crystals by con- 
trolling the detuning of a spin-dependent force in the 
region of closely-spaced transverse motional modes, and 



this approach should be scalable to much larger numbers 
of ions. One future direction is to add an effective mag- 
netic field (e.g., BJ^i'^y^) transverse to the Ising cou- 
plings, and characterize nontrivial phases of this system 
for larger numbers of ions [3, EBl • 
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